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Abstract 

We formulate the massive overlap fermions on anisotropie lattiees. We find that 
the dispersion relation for the overlap fermion resembles the continuum form in 
the low-momentum region once the bare parameters are properly tuned. The quark 
self-energy and the quark field renormalization constants are calculated to one-loop 
in bare lattice perturbation theory. Using tadpole improved perturbation theory, 
we establish the relation between the well-tuned parameters and the quark mass 
once the bare coupling is given. We argue that massive domain wall quarks might 
be helpful in lattice QCD studies on heavy-light hadron spectroscopy. 
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1 Introduction 



It is known that chiral symmetry plays an essential role in low-energy Quantum 
Chromodynamics (QCD), which is believed to be the underlying fundamental 
theory for strong interactions. Another predominant feature for low-energy 
QCD is its non-perturbative nature. Lattice QCD provides a genuine non- 
perturbative theoretical framework for the study of low-energy QCD from 
first principles. Chirality is associated with massiessness of fermions. Incorpo- 
rating massiess fermions have always been a great challenge in lattice studies 
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due to both theoretical and numerical difficulties. In recent years, consider- 
able progress has been made in understanding chiral symmetry on the lattice. 
Domain wall fermions [Tll^ and the overlap fermions piHI3]|5H7] have emerged 
as two new candidates in the formulation of lattice fermions which have much 
better chiral properties than the conventional Wilson or staggered fermions. 
Since chiral symmetry is so cruciai to the theory of QCD, it is therefore de- 
sirable to use these new fermions if possible. 

On the other hand, anisotropie lattices have been used extensively on heavy 
hadronic states and they proved to be extremely helpful in various appli- 
cations. These include: glueball spectrum calculations [8|9][T0] . charmonium 
spectrum calculations [TT][T2] . charmed meson and charmed baryon calcula- 
tions [rSfT^ and hadron-hadron scattering calculations [T5|16II17|[T8] . Note 
that many of the above mentioned studies involve light quarks for which chiral 
symmetry is essential. It is therefore desirable to use either domain wall or 
overlap fermions which have better chiral properties. Indeed, much of the sys- 
tematic uncertainties in these studies originates from chiral extrapolations. It 
is therefore tempting to study the new lattice fermions (domain wall fermions 
or overlap fermions) on four- dimensionai anisotropie lattices. In a previous 
study, we have formulated domain wall fermion on anisotropie lattices using 
bare perturbation theory to one-loop order [19]. In this paper, we will perform 
a similar study for the overlap fermions. These studies provide us with some 
guiding Information on the tuning of the parameters in the corresponding 
fermion action which is necessary in realistic Monte Carlo simulations. 

In this paper, we study massive overlap fermions on anisotropie four-dimensional 
lattices. For the gauge action, we adopt the tadpole improved gauge ac- 
tions [20][8] which have been used in various lattice calculations. The fermion 
action on anisotropie lattices generally contains more parameters than its 
isotropie counterparts. These parameters have to be tuned properly in or- 
der to yield a correct continuum limit. We will first address this issue in the 
case of free overlap fermions on anisotropie lattices. It is found that, in order 
to restore the normal relativistic dispersion relation for the quark, parameters 
of the fermion action have to be tuned accordingly. Then, we compute the 
quark propagator in lattice perturbation theory to one-loop. Quark field and 
quark mass renormalization constants are obtained for various values of the 
bare parameters. This perturbative calculation serves as a guidance for further 
non-perturbative Monte Carlo simulations. Similar perturbative calculations 
have been performed in the case of isotropie lattice [2T]|22|23||2^ . Our calcula- 
tion is an extension of these to the case of anisotropie lattice. The use of both 
anisotropie lattices and overlap fermions to treat relativistic heavy quarks on 
the lattice was first advocated in Ref. [25] where the authors have considered 
the dispersion relation of quarks in a quenched numerical study. 

This paper is organized as foUows. In section 2, overlap fermion action on 
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anisotropie lattices is given. In section 3, the free overlap fermion propagato! 
on anisotropie lattiees is presented and the dispersion relation of the free over- 
lap fermion is studied. It is found that hopping parameters of the fermion ac- 
tion have to be tuned properly, according to the vaine of the bare quark mass, 
so as to have the correct continuum limit for the massive quark. In section 4, 
the calculations of fermion self-energy to one-loop is presented. Numerical re- 
sults for the renormalization factors for the quark field and the current quark 
mass are listed. In section 5, we will conclude with some remarks and outlook. 
The one-gluon and two-gluon vertex functions are listed in the appendix. 



2 The overlap fermions on anisotropie lattices 



The fermion action for the exact massless overlap quark has the foUowing 
form PGIIT] : 

SF = Y^i,{x)D,yÌj{y), (1) 

x,y 

with the (massless) overlap Dirac operator given by: 

(2) 

The operator X appearing in the above equation is the Wilson-Dirac operator 
on anisotropie lattices which we write as: 
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(3) 



where represents Wilson parameters introduced to remove the doublers 
from low-energy spectrum. In the case of anisotropie lattice, we will also use 
the convention: tq = r^, = r^. The parameter Mq is usually taken to be 
negative so as to restore chiral properties of the quark. Since one can scale 
the operator X by a Constant factor without changing fermion matrix one 
might just fix the parameter to some Constant, say We will also set the 
temporal Wilson parameter = 1. It is well-known that the massless overlap 
Dirac operator Dira = 0) satisfies the Ginsparg- Wilson relation [26] which 
can be viewed as a generalization of the continuum chiral symmetry to the 
lattice [27f28] . The Dirac operator also has a massless mode which has definite 
chirality. 

For the massive quarks, one has to modify the massless overlap Dirac operator 
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given in ([2]) to [5]: 



It is known that the chiral mode then acquires a mass that is proportional to 
the parameter m for small values of m. With these conventions, the lattice 
action describing a massive domain wall fermion depends on 4 bare parame- 
ters: the Wilson mass parameter Mg, the temporal hopping parameter Kt, the 
spatial Wilson parameter and the bare quark mass parameter m. 

In perturbation theory, it is convenient to study these matrices in Fourier 
space. After Fourier transformation of the matrix X: 

the quantity X{q,p) is expanded into power series in the bare coupling go: 

X{q,p) = X,{p){2TiY5p{q-p) + X,{q,p) + X^iq^p) + 0{gl), (6) 

where the functions Xi are of the order of {g^Y respectively. Higher order con- 
tributions are not needed in an one-loop calculation. The explicit expressions 
are found to be: 



Xq{p) = + e '^f^'^À^ ~ cosp^) + M( 



0, 



^i(?.p) = E / ^05(2>r)''ip(<, -p-k)y g<,Aj(k)T^V,^ 



A,B,fj.,u 



. (27r)4 (27r)^ 
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where p = J2fi i^iilfi sinp^ and the functions Vi^ and are given by: 



Vi^,{p + -) = ft^ COS p + - + sin p + - 



V2JP + 



Eh 



/ 

Eh 



COS p 



Eh 
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3 Dispersion relation for the free overlap fermions 



In this section, we briefly summarize the resuhs for the free overlap fermions on 
anisotropie lattices. These results can be obtained from their counterparts for 
the isotropie lattices which can be found from the literature, see for example 
Ref . [5] • In the free case, we denote the overlap Dirac operator given in Eq. (jl]) 
by Dq. The inverse of Dq in momentum space is found to be: 



b{p) 



l(l-f)xt(p) + 


(l + f ) u{p) 


(9a) 


2(l + f)'^W + ( 


1 - f ) b{p) ' 


^K^r^(l - cosp^) 




(9b) 









2 



p^+ [Y.K^r^{l-cosp^) + Mo\ >0, (9c) 



where p"^ = J2fi ''^fi sin^p^. 

The dispersion relation for the physical particle (in this free quark) is 

obtained by inspecting the time dependence of the propagator in coordinate 
space: 

Go(p,t)= r ^D^\po,p)e'^°' . (10) 

Fot large temporal separation t, the above propagator will behave like e~*^p 
with Ep being the energy of the particle with three-momentum p. Using the 
variable z = 6*^°, the above propagator can also be expressed as: 

Go{p,t)= i^z'~'D-\z,p). (11) 
J Im 

where the contour integrai is along the unit circle counter-clockwise. If D~^{z^ p) 
were an entire function of z, the above propagator receives contributions from 
the poles of D~^{z) within the unit circle. The case for the overlap is a bit 
more complicated due to the branch cut of uj{p). However, it can be shown 
that the leading contribution to the propagator GqÌI, p) at large t comes from 
the residue of the pole closest to z = 1 in the complex z piane. Thus, setting 
the four-momentum: 

Pm = (^P'P)' (12) 
and demanding that it is the pole of the propagator Dq^{p), the above equa- 
tions yield the relation: 

_ (^i + Hpj'p^^rn'bK (13) 
We are only interested in the dispersion relation when the lattice three-momentum 
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IpI ^ 1. We thus obtain: 

where Eq = mg will be identified as the pole mass of the quark and Mkm is 
the so-called kinetic mass of the quark. After some calculations, the equation 
satisfìed by Eq = mq is found to be: 

A cosh^ £"0 + S cosh £^0 + = , (15) 



where the coefficients are given by: 



2\ 2 

2 / 1 I 1 „2_2 



A = K^m 



4 

,2 



B = iKtm' {Kt + Mo) , (16) 



,2' 
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It is seen that Eq = mg is independent of the parameter r^. It is also easy to 
verify that for small values of m, the solution Eq = mg obtained from Eq. f|T5l) 
is proportional to m, as expected. We can also find out the kinetic mass term 
with the result: 

1 _ dEp 

In order to have the usuai energy-momentum dispersion relation for the quark, 
one imposes the condition: Mkin = Eq = mq which yields the relationship 
between x ^"^^ '^t as (restoring the lattice spacing explicitly): 

^2 ^ (A cosh £0 + B/2) (sinh Eq/Eq) ^^^^ 
(1 + + mV, (Mo + Kt-Kt cosh £0) ' 

where the value of Eq = niq is to be obtained from Eq. f|T5|) . To summarize, 
for a given set of parameters Mq, m and r<., one has to solve both Eq. f|T5l) 
and Eq. fITH]) to obtain the values of Kt and mq. 

In Fig. [H we have shown the pole mass of the quark as a function of the 
propagator mass parameter: mp = m|Mo| for a given set of other parameters. 
Three curves in the plot correspond to x = I7 3, 5 respectively. The value of Kt 
is obtained from Eq. (fTSl) . It is seen that the pole mass of the quark increases 
with rrip linearly for small values of mp. For larger values of bare quark mass, 
some non-linearity sets in. 



r>=n 



(1 + ^) + mV, (Mo + Kt-Kt cosh Eq) 
{2A cosh £0 + 5) sinh £0 



(17) 
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Fig. 1. The pole mass of the free domain wall quark, mg, measured in l/os unit as a 
function of the propagato! bare mass parameter mp = m|Mo| for three values of the 
anisotropy parameter x- Three curves correspond to x = 1, 3 and 5, respectively. 
We have set Mq = —0.5 and = 1 in this plot. 

In Fig. [21 for a given value of anisotropy Xì w6 have plotted the appropriate 
vaine of Kt as a fnnction of the propagator mass parameter mp. Since mp 
is very dose to the pole mass mg as Fig. [1] indicates, this fignre can also be 
viewed as the dependence of Kt on mg for a given anisotropy x- K is seen that, 
in the massless limit, i.e. mg — ^ 0, one recovers the naive relation: x = '^t- 
However for massive quarks, this relation is distorted with increasing values of 
the quark mass. The deviation from its massless limit can be as large as 15% 
for values of mp ~ 0.5. In quenched Monte Carlo simulations, the anisotropy x 
is fixed by the pure gauge sector. Therefore, this figure can be utilized to tune 
the hopping parameter Kt accordingly for a given value of the quark mass and 
the anisotropy. In a non-perturbative Monte Carlo simulation, this tuning 
process can be performed by demanding physical particles, e.g. the pions, 
have the correct dispersion relation for small three-momenta. We would like 
to emphasize that this tuning process is cruciai in the simulation of massive 
quarks since, without it, the hadrons will not have the correct continuum limit. 

In both Fig. [T] and Fig. [2] we have chosen the Wilson parameters to be unity. 
In principle, other values are also permitted as long as the doublers are well 
separated from the physical modes. It is seen from Eq. (19bp that the separation 
between the doublers and the normal physical quark is characterized by the 
Wilson parameter r^. Therefore, to decouple the doublers it is better not to 
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Fig. 2. With the sanie set of parameters as Fig. [H the value of t^tlx is shown as a 
function of the propagator mass parameter mp for x = 1) 3 and 5. It is seen that 
in the massless hmit, the expected relation: x = is recovered for ali x- 

take too small values for r^. In the remaining part of this paper, we will take 
the conventional choice: = = 1. 



4 The wave function and quark mass renormalization to one-loop 



In this section, we will compute the quark self-energy to one-loop using bare 
perturbation theory. Two diagrams contribute at this level: the tadpole di- 
agram and the half-circle diagram. We will first list the necessary Feynman 
rules and the calculation of these two diagrams will be dealt with afterwards. 



^.i Feynman rules 



First of ali, one needs the free propagator of the lattice gauge fields [29] . After 
performing Fourier transformation of the gauge fields, the quadratic part of 
the gauge action has the standard form in momentum space: 



si'^lA,] = ^ E £ 7^ (MI) MAO M-i)) , (19) 
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where 



Moo=x 



1 

X 
1 

X 



r 



a, 



a, 



a, 



■^0 + E??oj 



7 2 ,2/2 , i2 



Moj=M,o = X 



— IqIj — lolj qoj 



L"9 



(20) 

(21) 

(22) 
(23) 



with lattice momentum defined as: = 2 sin(y). The quantities appearing 
in the above equations are given by: 



Using these notations, the free gluon propagato! is expressed as: 



(24) 



[V 

The exphcit expressions for f^^ and fo maybe found in the hterature [29]. For 
simphcity, in the foUowing calculation we choose the gauge in which Ug = 1. 

The vertex functions for the interaction between the quark and the gluon 
fields can also be obtained with the help of Eq. ([7]) and Eq. ([8]). The explicit 
expressions are given in Eq. (lA.ip and Eq. (1A.2P in the appendix. 



r{ip.qp.:x) 

foilo.Qou.X) . 



(25) 



4-2 The half-cirde diagram 



For the half-circle diagram, the contribution to the quark self-energy can be 
written as: 



rl^h 1 — — 

^half-circlC|'^\ _ J2n 1 \ 2 



X { -7 ( n (p _ k) , (26) 



where 
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1 



2 ; uj{p) 2 J uj{k) 



and Dfj^fj^ is the gluon propagato! given in Eq. fl25l) . 



We will be interested in the small p behavior of the self-energy. Therefore, 
one can expand the integrand in Eq. fl26l) around p = 0. Note that uj{p) is 
an even function of p. Therefore, when expanded around p = 0, one has: 
u;{p) = uj{0) + 0{p'^). However, both D^^{p — k) and crj^^{p, k) contains terms 
that are hnear in p. When expanding the function D^^{p — k), the leading 
terni is an even function of k, the terni proportional to p is an odd function of 
k. Therefore, when multiphed with expansion of the function a^^\p,k), only 
the terms that are even functions of k will give non-zero contributions after 
the integration over k is performed. Por example, if we take Z)^^(/c), then it 
suffice to keep only even functions of k in the expansion of (T^^\p, k): 



a«(0,A;) 




(28) 



Here the symbol "~" means that the odd functions of k are dropped from the 
expression. Similarly, if the terms linear in p from D^^{p — k) are taken, which 
is an odd function of /c, then only terms that are odd functions of k in the 
expansion of a^^^ (p, k) will contribute to the integrai. 



Then, we calculate the quantity da^^\p, k)/dp^. First, we notice that 



dp^ 



p=0 



dp^ 



(29) 



p=0 



We define: da^\p^ k)/dp^ = in^'j^a^^^''' and the explicit calculation shows: 
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p=0 



Ut 



1^ — ^1/1" 



z . (k\ Xo{k) ~ 



Lj(A;) 



vi I ^ 



X 



X 



2 

m 



lu 



Xoik) . (k\ i 



a;(A;) 



(30) 



As we clarifìed before, we only need to keep odd functions of k in this quantity. 
Thus we get: 



da^^){p,k) 



dp^ 



p=0 



m 



2 



'-{^,, V,,XlVi,}5,u + Vlx,vl}5,. 



Il 



m 



+ [1 + f j \ -^1,^ [Vi,V,, + vivi) 5,, 

vlXoVl}5,. + i^.^h iylv,. + ì/i.vit) 



Mo u 



(31) 



Here the symbol "~" means that the even functions of k are dropped from 
the expression. Using the expressions for various quantities, we have: 



ViyViy ~ -kI cos ky vivi ~ -kI cos ky 
VIXqVI ~ kI {k^ sìv? ky-h cos kl) 

ViuX^Viu ~ (^K^ sin^ /^i/ - è cos /e,,) 
Finally, we calculate the terni dD^^{p — k)/dp^\p=o: 



(32) 
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dDy^{p - k) 



p=0 



1 



JD 



-Cu -CD £vu 
J u , J u J 



(33) 



where we bave used the foUowing notations: 



df^ip-k) 



k^dp^ 



f^ik) 



dfnip-k) 



p=0 



k^^dpf. 



p=0 



(34) 



As mentioned above, here the even functions of k in ^"^^■'(O, k) will vanish and 
the odd ones will survive. These terms are of the foUowing forni: 



k^a^^^ (0, k) ~ iKf.'y^k^ sin /c^ <( 2 ( 1 + y ) K^ujóf,^ 



fyi\ l 2 „;„2 ky I 2 2 ky 



+2 ( 1 + - j sin^ -^ + KÌ cos' ^ {25^y - 1] 



m 



a; V 2 



-bn^S^y + ni- 2kI cos^ ^5^^ 



(35) 



Collecting ali the relevant terms, we may write the contribution from the 
half-circle diagram as: 



-ihalf — circlc 



Va* •> 



(36) 



Introducing the foUowing notations: 



da(J\p,k) 



dPt. 



tK,^,a^^>{k) k,a^J\0, k) = zK,^,àl^>{k) , (37) 



p=0 



the quantity J*-^-* and M^^^ are expressed in terms of the foUowing integrals: 
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(2ir)4 ^ - MoY 
A\ , f^f U,.,W. (39) 

where ali integrations are taken in the first Brillouin zone. 
4-3 The tadpole diagmm 



Using the Feynman rules given in the previous section, the contribution from 
the tadpole diagram can be written as: 



Stadpoie(p) = \g'C, £ ^ E (l - y) ^f\P^ k)D,,ik) , (40) 
where the function crj^^p, k) reads: 



2uj{p){ '^^"^ uj{p) uj{p) j uj{p){uj{p)+uj{p + k)y 

X ixoipWl (p + ^] V,, (p+^]+ VI, (p + xUp + A:) VI, f |> + ^ 



+V^,[p + ^)vl(p+^]Xo{p) 



'^§^MPWi(p.^^ 

(41) 



For small lattice momenta p, the function crj^\p, k) is expanded. It contains a 
terni at vanishing moment um: 



(0, A;) = (ylV,, + + -nl^Xo^ij , (42) 
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and a terni that is linear in p. To evaluate this terni, we need da^\p, k)/dp^ 



daf\p,k) 



dPn 



p=0 



' Mo 



\ . 4:kI cos kf, + 4k^6 
^ Mo / '"^^'^ {cu - Mo^Mou^ 



sin k,. 



X I -hK^ó^u + 1^1- cos^ y^^^ 
1 



Mo{c.'(A;) - Mo}2 



- z{7^, Vi^XoVi^j^^^ + Fi 



dPt. 



p=0 



2fi;2 cos k^ + 2k^6 . 

— ^ sm k^ 



X Mo + '2^1 cos^ ^^/.;. 

1 2 



Mo 



+ Mont 



dX o{p + k) I 
dp^ 



lp=0^1i/ 



(43) 



Note that we have: 

Vi. 
Vi 



dXl{p^ 


-k) 


Vi. ~ 

p=0 


dp^ 




dXoÌp + 


k) 


vl = 

p=0 


dpf. 



Vii, ~ -ii^tilfi cos /c^K^ + mj^7^(l + cos k^)6^ 



dxUp + k) 

dpf. 



\p=o 



Vi. 



(44) 



Therefore, if we introduce the foUowing denotation 
da(^\p,k) 



dp^ 



(45) 



p=0 



we can parameterize the contribution from the tadpole diagram as: 

E-'^-^^(p)=/(^.E7.;5,jf +M(^)) , 
where the quantities /'•^^ and M*^^^ are given by: 



(46) 
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= 7^^ E fi - T ) -?''ik)DUk) , (47) 



m 

> 2^"y_. (27r)4^ 2". 



Adding the contributions from the half circle and the tadpole diagrams to- 
gether, one can verify that M^^^ + M^'^^ is proportional to the quark mass: 



2 + 



^mm(^) • (49) 



This is the expected result since it is known that quark mass is only multi- 
phcative renormahzed in QCD due to chiral symmetry. In what foUows, we 
will define: 

M(^) = , 50 

mp 

which is a finite quantity in the chiral hmit. 



4-4 The massless case 



li the quark mass is exactly zero, then the one-loop contribution discussed 
above contains infra-red divergences. However, the same infra-red divergence 
also exists in the continuum. Therefore, one can subtract an appropriate infra- 
red divergent part from the one-loop contribution: 



,(1) _rii)f _a^ cp r d'k ^x{x'Y^%o{^' - {xX + J:,kjl\ , 

Sfinite,/. -V l"^-Uj / _ ^^^^4 (X^'^O +Ei A;|)3 




'S-. = T^|Mvr^) + l-Ml5')l (52) 



Mi:' = -T-f^T H^') + 1 - HP') (54) 
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The subtracted contributions listed above are ali infra-red finite. Therefore 
they can be evaluated numerically once the bare parameters are fixed. The 
infra-red divergent part can be obtained analytically and depends on the exter- 
nal momentum (or a scale parameter when external momentum is vanishing) 
as usuai. 



4.5 Numerical resulta 
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Fig. 3. The values of \m) as a function of mp are shown for x = 1, 3 and 5. Here 
/j^^ corresponds to ^ = and corresponds to /U = 1, 2, 3. Other bare parameters 
are: Kg = 1, Mq = —0.5. 

As we pointed out above, if the quark mass is non- zero, quantities lj^\m), 
lj^\m) and M^^\m) can be calculated directly using numerical integration 
once the bare parameters are given. In Fig. [3] and Fig. HI we have shown the 
values of ij^^ and ij^^ as functions of the propagator mass parameter mp. Here 
the subscript "t" corresponds to /i = and the subscript "s" corresponds to 
fi = 1,2,3. In Fig. we have shown the values of M^^^ as a function of the 
propagator mass parameter rrip. 

In the massless case, the subtracted part of the loop integrals can be computed 
numerically foUowing Eq. (!5T]) and Eq. (!53l) . Since the main purpose of this 
paper is to address massive overlap quarks, we will not list the numerical 
results for the massless case. 

Finally, our calculation can be easily translated into its tadpole-improved ver- 
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Fig. 4. The values of (m) as a function of mp are shown for x = 1, 3 and 5. Here 
corresponds to ^ = and 1^^^ corresponds to /U = 1, 2, 3. Other bare parameters 
are: = 1, Mq = —0.5. 

Sion foUowing standard steps. For example, if we use the mean-field estimate 
(tree-level tadpole improved theory), Eq. (l9bll is modified to: 

b{p) = ^i^f,{l - u^icospf,) + Mq , (55) 

where represents the mean-field vaine for the gauge field t/^. Note that for 
small lattice momenta, this amonnts to a shift in parameter Mq: 

Mo = Mo + ^/€^(l-n^) . (56) 

Similarly, the tuning of the hopping parameter Kt can also be discussed within 
the mean-field approximation. Note that these might be a rather good estimate 
for the correct vaine for the parameters in future Monte Carlo simulations. 



5 Discussions and conclusions 



In this paper, we have studied massive overlap fermions on anisotropie lattices. 
We argue that this setup can be useful in many lattice QCD studies. It is shown 
that, in order to restore the usuai dispersion relation for the quark at small 
three-momenta, hopping parameter has to be tuned according to the quark 
mass values. Quark propagator is calculated using bare perturbation theory to 
one-loop order. We find the wave-function and mass renormalization constants 
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Fig. 5. The values of M^^\m) as a function of mp are shown for x = 1) 3 and 5. 
Other bare parameters are: = 1, Mq = —0.5. 



at various values of the bare parameters. These results serve as a guidance for 
the tuning of the parameters in real Monte Carlo simulations. 



A Vertex functions 



Here we list the vertex functions used in the main text. The isotropie lattice 
counterparts for the vertex functions can be found in the literature [2T|22f23|2^ 
and the version for anisotropie lattice can be obtained accordingly after obvi- 
ous modifìcations. The interaction vertex for the quark field, anti-quark field 
and one gluon field reads: 



uj{q)+uj{p) [ \ 



The the two-gluon interaction vertex is much more complicated. It is found 
to be: 
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uj{q)+uj{p) 



V2JP 



{uj{q) + uj{p)}{uj{p) + uj{p + ki)}{uj{p + ki) + uj{q)} 
y<[xo{q)Vl[p+k, + ^y,^(^p+fj 

+ Vl, (^P+k2 + X^,{p + k2)V,,{p + |) 

+Vi, (^p+k2 + y ^ vi (^p + Xoip) 

uj{q) + uj{p) + uj{p + A;2) ^ 
uj{q)uj{p)uj{p + k2) 

M<l)Vl [T> + k2 + Xo(p + k2)Vl [p + Xo(j))| 

1 

+ u;(|})}{u;(p) + uj{p + A;i)}{u;(p + ki) + 
iv + h + ^] + (p + y ] 



k. 



Vi P 



2 ) '^^y" ' 2 
uj{q) + + u;(p + ki) 



u{q)uj{p)uj{p + ki 



X 



(A.2) 
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